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We show that layered quenched randomness in planar magnets leads to an unusual intermediate 
phase between the conventional ferromagnetic low-temperature and paramagnetic high-temperature 
phases. In this intermediate phase, which is part of the Griffiths region, the spin-wave stiffness per- 
pendicular to the random layers displays anomalous scaling behavior, with a continuously variable 
anomalous exponent, while the magnetization and the stiffness parallel to the layers both remain fi- 
nite. Analogous results hold for superfluids and superconductors. We study the two phase transitions 
into the anomalous elastic phase, and we discuss the universality of these results, and implications 
of finite sample size as well as possible experiments. 



PACS numbers: 75.10.Nr, 75.40.Cx, 74.40.-n, 67.85.Hj 
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The macroscopic behavior of many-particle systems is 
often sensitive to quenched disorder. For example, at 
zero-temperature quantum phase transitions, the inter- 
play of quantum and disorder fluctuations gives rise to 
exotic phenomena, such as quantum Griffiths singulari- 
ties [1-3], infinite-randomness critical points [4, 5], and 
smeared transitions [6, 7]. The main reason for these 
strong effects of disorder is the presence of perfect disor- 
der correlations in the imaginary-time dimension, which 
becomes infinitely extended at zero temperature. Thus, 
one is effectively dealing with infinitely large impurities. 

This suggests that strong disorder effects should also 
occur at classical (thermal) phase transitions, if the dis- 
order is perfectly correlated in one or more space dimen- 
sions. For example, the McCoy- Wu model [8, 9], a two- 
dimensional (2D) Ising model in which the disorder is 
perfectly correlated in one dimension, shows an exotic 
transition, characterized by a smooth specific heat but 
an infinite susceptibility over a range of temperatures. 
By using a strong-disorder renormalization group, Fisher 
[4, 5] showed that the critical point is of the infinite- 
randomness kind, and is accompanied by power-law Grif- 
fiths singularities. Similar behavior was found in Heisen- 
berg magnets having 2D disorder correlations [10]. 

In this Letter, we study thermal phase transitions ex- 
hibited by randomly layered 3D superfluids, supercon- 
ductors, and planar magnets, as sketched in Fig. 1. All 
these systems are characterized by two-component or- 
der parameters of U(l) or, equivalently, 0(2) symmetry 
(representing the condensate wave function, Cooper pair 
amplitude, and magnetization, respectively). 

Couched in terms of the planar ferromagnet, our re- 
sults can be summarized as follows: The interplay of the 
layered randomness and the Kosterlitz-Thouless (KT) 
[11] transitions in strongly coupled multilayers (slabs) 
leads to an anomalously elastic intermediate phase as 
part of the Griffiths region associated with the phase 
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FIG. 1: (Color online) Schematic behavior of the magnetiza- 
tion m and the stiffnesses p„n and p St ± vs. temperature T 
for a bounded disorder distribution. SD and SO denote the 
conventional strongly disordered and ordered phases, respec- 
tively. The Griffiths region (bounded by T u and Ti) consists 
of the "non-anomalous" (G) and the anomalous (AG) Grif- 
fiths phases. For an unbounded distribution, T u — > oo. Inset: 
Randomly layered magnet or superconductor: layers of two 
distinct materials are deposited in a random sequence. 



transition. In this anomalous Griffiths phase, the mag- 
netization m and the spin- wave stiffness p Sj y parallel to 
the layers are both nonzero (as in a conventional ferro- 
magnet). However, the stiffness p s .± perpendicular to 
the layers vanishes, and the elastic free energy exhibits 
anomalous scaling behavior. Specifically, free energy AF 
due to twisted boundary conditions (BCs) in the z direc- 
tion, forcing the spins on the top face in Fig. 1 to make 
an angle of O with those on the bottom face, varies as 
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with system size L±. Here, z(T) is a temperature depen- 
dent dynamical exponent that varies continuously from 
z = oo at the upper Griffiths temperature T u (i.e., the 
boundary between the Griffiths region and the conven- 
tional paramagnet) to z = 1 at the temperature T s where 
a nonzero perpendicular stiffness first appears. 
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While anomalous elasticity of this type occurs in some 
disordered systems possessing uncorrelated disorder (e.g., 
liquid crystals in aerogel [12]), in those cases it is charac- 
terized by universal values of the anomalous exponent z: 
the non-universality and strong temperature dependence 
of z that we find here are, to the best of our knowledge, 
unique to systems having correlated disorder. 

We also find unusual behavior at the two transition 
temperatures T u and T s . The magnetization m is nonzero 
for all T < T u and shows a double-exponential tail to- 
wards the nonmagnetic phase. Close to T u , it takes the 
asymptotic form 

ln(m) ~ - cxp[a(T u - T)~ v ] , (T -> T u -) (2) 

where v k, 0.6717 [13] is the correlation-length critical 
exponent of a clean 3D planar (XY) magnet and a is a 
nonuniversal constant. If an external magnetic field H is 
applied at temperatures T <T Ul the magnetization van- 
ishes with decreasing field more slowly than any power, 

ln(m) ~ - VI ln(ff) | (T u -T)-» , (H -> 0) . (3) 

This relation applies for magnetizations larger than the 
double-exponentially small value given in (2). 

The parallel spin-wave stiffness (corresponding to a 
twist of the BCs in the x or y direction) p s u is nonzero 
for all T <T U and shows an exponential tail of the form 

In(p a ,||)~(T u -T)- v , (T^T U -). (4) 

In contrast, the perpendicular stiffness p Si ± (correspond- 
ing to a twist of the BCs in the z direction) vanishes as 
T s is approached from below via 

/VL ~ (T a - T) , (T^T S 7). (5) 

In the remainder of this Letter, we sketch the deriva- 
tion of these results, compute finite-size effects, and dis- 
cuss possible experimental realizations. For definitcness, 
we focus on the classical planar ferromagnet (i.e., the XY 
model) on a cubic lattice. The Hamiltonian is given by 

H = — ^ jj (S r -S r+ x+S r -S r+ y)— ^ J^r S r -S r+ z. (6) 

r r 

Here, S r is a two-component unit vector at lattice site 
r, and x, y, and z are the unit vectors in the coordinate 
directions. The exchange interactions within the layers, 
jj, and between the layers, J^, are both positive and 
independent random functions of the perpendicular co- 
ordinate z. For simplicity, we take all Jjr = J 1 - and draw 
the jj from a binary distribution (J u > J{) 

P(JH) = (1 - c) 5(jH - J u ) + c6(jW - Ji) , (7) 

where c is the concentration of "weak" layers. 

Let us discuss the planar magnet (6) qualitatively. At 
sufficiently high temperatures, the system is in a con- 
ventional (i.e., strongly disordered) paramagnetic phase. 



Below the upper Griffiths temperature T u , which is de- 
fined as the transition temperature of a clean system hav- 
ing Jz = J u , rare thick slabs ("rare regions") of strong 
(Jz = J u ) layers show local magnetic order, while the 
bulk is nonmagnetic. Although individual such slabs 
are prevented from developing true long-range order [14], 
they can undergo KT transitions [11]. The unusual be- 
havior, eqs. (1) to (5), is caused by the interplay between 
the randomness and the KT physics of the rare regions. 
Ultimately, below the lower Griffiths temperature Tj (the 
transition temperature of a clean system having jj = J; ) , 
all layers order magnetically, and the system exhibits the 
conventional (i.e., strongly ordered) ferromagnetic phase. 

We now use optimal fluctuation theory (i.e., Lifshitz- 
tail arguments [15]) to derive the thermodynamics in the 
Griffiths region. The probability w(L RR ) f° r finding a 
rare region of L RR consecutive strong layers reads 

w(L RR ) ~ (1 - c) Lrr = e~' cLRR , (8) 

with c = — ln(l — c). Each individual such slab is 
equivalent to a 2D XY model, and thus undergoes a 
KT transition at some thickness-dependent temperature 
Tkt(Lrr). Finite-size scaling yields T u - T KT (L RR ) ~ 
Lrr ■ This result defines a cutoff length L C (T) <~ 
(T u — T)~ u . At any temperature T < T u , all rare regions 
of thickness L RR < L C (T) are (locally) in the disordered 
phase, while those having L RR > L C (T) are in the quasi 
long-range ordered KT phase. 

Let us first consider a single rare region. According to 
KT theory [11], the spatial correlation function C(x) in 
the KT phase falls off as a power of the distance |x|: 

C(x) ~ |xp (|x| -+ oo) . (9) 

The exponent 77 is related to the renormalized (parallel) 
spin- wave stiffness p s , RR of the slab via n — T / (2irp StRR ). 
It takes the value 1/4 at the KT transition and is 
inversely proportional to L RR for very thick rare re- 
gions. We thus model the thickness-dependence of 77 via 
77 = jL c (T)/L RR , which correctly describes both lim- 
its. The power-law correlations also lead to a nonlinear 
magnetization- vs. -field curve within the KT phase, 

m - iTJ/( 4 -") , (10) 

which implies an infinite magnetic susceptibility. 

We now combine the single-slab results, (9) and (10), 
with the size distribution (8). We start with the response 
to an external magnetic field H. Neglecting interactions 
between the rare regions for the moment, we write the 
rare-region contribution to the magnetization as 

/>oo 

m~ / dL RR w(L RR ) H^ L ^y^ L ^ . (11) 
Jl c (t) 

For small fields, this integral can be evaluated using the 
saddle-point method. This yields a saddle-point value 
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L 2 p = | ln(H)\L c (T) / (16c) , implying that the response at 
small fields is dominated by thick rare regions. Inserting 
L sp into (11) immediately gives (3). This highly singular 
result breaks down for H > H x ~ exp[— 16cL c (T)], at 
which the saddle point reaches L C (T). For H > H x , the 
response is of the power-law type, to ~ if 1 / 15 , until the 
regular, linear-response part takes over at an even larger 
field H rcg ~ exp[-(15/14)cL c (T)] [16]. 

The parallel spin- wave stiffness p s u can be found anal- 
ogously. The free energy cost due to a twist of the BCs in 
either the x or y direction is simply a sum over all slabs 
in the KT phase. Each slab has the same twisted BCs, 
thus, the total parallel stiffness is given by 

r-OO 

Ps,\\ ~ / dL RR w(L RR ) Ps,rr{Lrr) . (12) 

JL C (T) 

This integral is dominated by the contribution near the 
lower limit, where p s , RR is approximately constant and 
equal to 2T/tt. To leading exponential accuracy, we thus 
obtain p s n ~ exp[— cL c (T)], which leads to (4). 

To discuss the perpendicular stiffness p s ,±, we apply 
twisted BCs in the z direction. The resulting local twists 
occur mostly in the disordered bulk between the rare 
regions; due to the randomness, they are not uniform 
but vary from layer to layer. As the spatial positions 
of the rare regions are completely random, the distribu- 
tion of their nearest-neighbor distances R is Poissonian, 
P(R) = R K TCxp(-R/R KT ), where R KT ~ exp[£L c (T)] 
is the typical separation. The effective coupling be- 
tween neighboring rare regions falls off exponentially, 
J^g(R) ~ cxp(— R/£o), where £o is the bulk correlation 
length. Combining this exponential form with P(R) gives 
a power-law distribution for the effective couplings, i.e., 

P(J^) ~ (J^i- 1 . (13) 

The Griffiths dynamical exponent z = Rkt/£,o takes the 
value oo at T u , and decreases with decreasing tempera- 
ture. Writing the free energy due to the twist of the BCs 
by an angle 9 as AF ~ J2 Z -4s ©*> with J2 Z = ©. 
and minimizing w.r.t. the Q z , we obtain [10] 

p B , ± ~ (l/Jei)" 1 (14) 

where (• • • ) denotes the average over the distribution 
(13). This average diverges for z > 1, implying p s _±_ = 
at temperatures just below T u . Upon lowering T fur- 
ther, the exponent z reaches the value 1 at a temperature 
T s < T u . For T <T S (i.e., z < 1), the average converges, 
yielding a nonzero stiffness. Close to z = 1, the average 
behaves as (1/J^) ~ 1/(1 — z) yielding (5). 

Finally, we turn to the spontaneous magnetization to. 
The reason that to > for all T < T u is the infinite 
susceptibility of those slabs that are in the KT phase. 
They align to one another via an infinitesimal coupling. 
In contrast, in the quantum Griffiths scenario, realized 



in the layered Hciscnberg magnet [10], the rare regions 
have a large but finite susceptibility. Aligning them re- 
quires a nonzero coupling, so that long-range order only 
appears at some critical temperature below T u . To esti- 
mate to, we combine the effective interaction with 
the KT scaling within the rare regions. Consider an 
area of linear size L (in the x and y directions) in one of 
the slabs. The typical magnetization (per site) of such 
a region can be calculated by integrating (9), yielding 
m(L) ~ L~ v / 2 . Now consider two such areas in neigh- 
boring rare regions. Their interaction can be estimated 
as J^(L) = J^L 2 m 2 (L) ~ L 2 ~ v exp(—R KT /£,o)- When 
this interaction becomes of order T, the areas align, 
and long-range order sets in. This happens at a length 
L = L X ~ [exp(iJ Jf T/&)] 1/(2 -" ) , yielding 

to ~ L x ^ 2 ~ exp[-(iW£o)??/(4 - . (15) 

Because of the exponential size distribution (8), the vast 
majority of rare regions in the KT phase are very close 
to the KT transition. Thus, to a good approxima- 
tion, we can set r\ = 1/4. Inserting this, along with 
Rkt ~ exp[cL c (T)], into (15) yields the final result (2). 
This calculation can be refined by taking into account 
the random distribution of rare-region separations, which 
only modifies the nonuniversal constants in (2) [16]. 

We now turn to the aspects of finite system size. The 
main effect of a finite perpendicular size L±, which is 
experimentally important because the number of layers 
in a real sample will often be small, is to limit the max- 
imum rare-region thickness L^ R X in the sample. Esti- 
mating Lft R K via the condition that a sample of size L± 
contains, on average, exactly one such rare region, i.e. 

l ± w ( l rr) ^ L we obtain l rr ~ H L ±)/c- 

We note that L^ff introduces an upper limit to the 

integral (11) for the m(H) curve. When the saddle-point 
value L sp is larger than L^ R K , which happens for fields 
H <H X with ln(H x ) - ln 2 (L ± )/[cL c (T)], the integral is 
dominated by the contribution near the upper limit. For 
very low fields, (3) gets thus replaced by a power law with 
a size-dependent exponent: to ~ u BcL ^( T )/ ln ( L ±) ^ w ^ n 
B a constant. The same mechanism also introduces an 
upper limit into the integral (12) for the parallel stiffness. 
As this integral is dominated by the lower limit, the finite 
size only matters when L C (T) > L RR X . Thus, the expo- 
nential tail (4) of gets cut off near the upper Griffiths 
temperature, for T u - T < \ln(L ± )/c\- 1 / v . Using (13), 
the minimum in a sample of size L± behaves as L^_ z . 
Inserting this into the elastic free energy expression given 
above (14) yields the anomalous elasticity scaling (1). 
As an example of the effects of a finite in-plane size 
we discuss the magnetic susceptibility. When L\\ is 
finite, the susceptibility of a single slab in the KT phase 
is no longer infinite. Its L\\ -dependence can be obtained 
from integrating (9) to an upper cutoff Lu, which yields 

Xrr(L\\) ~ L,, ~ v . Summing this over all rare regions, 
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and evaluating the integral in the saddle-point approxi- 
mation, gives a total susceptibility (per unit volume) of 

In summary, we have shown that the randomly lay- 
ered planar magnet features anomalous elasticity and un- 
usual thermodynamics in parts of the Griffiths phase. 
Although we have considered the binary disorder distri- 
bution (7), the functional forms of the results (1) to (5) 
remain valid for any bounded distribution, provided it 
does not vanish too rapidly at the upper bound. If the 
distribution is unbounded, the tails of magnetization and 
parallel stiffness would extend to T = oo, implying that 
the system is always in the magnetic phase [16]. 

Our theory describes the regime where the system con- 
sists of a few isolated rare regions in a disordered bulk; it 
becomes controlled for T — > T u . To describe the forma- 
tion of bulk order close to 7], the growths and merging of 
rare regions need to be included. Moreover, the charac- 
ter of the vortex unbinding transition changes for layers 
that are coupled to already ordered slabs [17, 18]. 

The results (1) to (5) have been formulated in terms 
of the planar ferromagnet. Nonetheless, they apply to all 
transitions having 0(2) or U(l) order parameters, if ex- 
pressed in terms of the appropriate variables. For layered 
superfluids and superconductors [25], the magnetization 
should be exchanged for the condensate wave function 
or the Cooper pair amplitude, respectively. In the same 
way, the spin-wave stiffness should be exchanged for the 
superfluid density, and the external field could possibly 
be realized via the proximity effect. 

Let us relate our theory to the classification of phase 
transitions with disorder based on the rare-region dimen- 
sionality <1rr [7, 19]: It states that the critical behavior 
is conventional if dun is smaller than the lower critical 
dimension d~ of the corresponding clean transition; if the 
rare regions order independently (i.e., if dun > d~), the 
transition is smeared. The marginal case, dan = d~ , 
usually leads to an infinite-randomness critical point. 
Based on these arguments, one might expect an infinite- 
randomness critical point in our system. However, the 
quasi long-range order that arises on rare regions in the 
KT phase actually leads to a hybrid between a smeared 
and a sharp transition. On the one hand, long-range or- 
der is present in the entire Griffiths phase (extending to 
T = oo for an unbounded disorder distribution), just as 
at the smeared transition of the randomly layered Ising 
model [20]. On the other hand, the long-range order is 
due to a collective effect (rather than individual freezing 
of rare regions), as in the randomly layered Heisenberg 
magnet [10], which has a sharp transition. 

Not only are our results of conceptual importance for 
the theory of phase transitions, but also they can be 
tested experimentally by producing layered nanostruc- 
tures of magnetic or superconducting materials. Mag- 
netic multilayers having systematic variations of T c from 
layer to layer have recently been produced [21], and our 



theory should describe random versions of such struc- 
tures (with XY spin symmetry). Moreover, using ultra- 
cold atomic gases, one should be able to completely en- 
gineer the appropriate many-particle Hamiltonian. We 
note that the Kostcrlitz-Thoulcss transition in a single 
slab of an 87 Rb gas has already been observed [22]. 
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0339147, DMR-0906566 and DMR-0906780, and by the 
Research Corporation. We also acknowledge the hospi- 
tality of the Aspen Center of Physics. 

After completion of this work, we learned of a study of 
the same issues by means of a numerical strong-disorder 
renormalization group [23]. Our phase transition sce- 
nario agrees with that of Ref. [23], and our asymptotic 
analytical results complement their numerical data. 
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